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Abstract 

We compute the next-to-leading order QCD corrections to the polarized (and unpolarized) 
cross sections for the production of a hadron accompanied by an opposite-side prompt photon. 
This process, being studied at RHIC, permits us to reconstruct partonic kinematics using exper- 
imentally measurable variables. We study the correlation between the reconstructed momentum 
fractions and the true partonic ones, which in the polarized case might allow us to reveal the 
spin-dependent gluon distribution with a higher precision. 
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I. INTRODUCTION 

It is well known that only a small fraction of the spin of the proton is carried by quarks. This 
information was gathered through the study of longitudinally polarized deep-inelastic scattering 
(DIS) processes. However DIS data is not enough to extract the full x-dependence of the polarized 
quark and gluon densities of the nucleon Furthermore, the x-shape of Ag seems to be hardly 
constrained by DIS data because it contributes in the leading order (LO) only through the Q 2 
dependence of g\. 

The precise extraction of Ag thus remains one of the most interesting challenges for spin physics 
experiments. Unlike DIS, processes involving the collision of polarized protons, like at the RHIC 
collider at BNL, introduce a direct gluonic contribution already at the lowest order in perturbation 
theory, resulting in a very useful tool to extract more information about Ag. 

A recent next-to-leading-order (NLO) global analysis, performed including all available inclusive 
and semi-inclusive DIS and RHIC data, suggests that the gluon polarization in the nucleon is 
rather small in the region 0.05 < x < 0.2 [l| Q|. Measured observables in hadronic collisions 
include single-pion production at center-of-mass energy yfs — 62 Q and 200 GeV [4|, and jet 
production at y/s = 200 GeV [5j. 

Both experimental collaborations at RHIC have started to study less inclusive observables that 
might turn out to be useful to pin down the shape of the polarized parton distributions with better 
precision. For example, the production of a charged hadron accompanied by a back-to-back jet 
or prompt photon allows experimentalists to use the last one as a trigger for hadron detection, 
reducing the bias in the selection of events [6j. Furthermore, from a theoretical point of view, 
counting with an extra particle (like a direct photon) in the final state provides a cleaner resource 
to access parton kinematics. The dominance of the partonic qg — > channel in the case of photon 
production supplies a direct signal on the gluon distribution in polarized collisions. 

Looking at a photon as a trigger has a relevance also in heavy ion collisions. When a hot and 
dense high-energy medium is formed, QCD particles interact strongly with this medium, causing 



them to lose energy as observed in Refs. Since the photon weakly interacts with that 

medium, this particle might travel through it without being strongly perturbed. For that rea- 
son, observing direct photons in coincidence with hadrons turns out to be useful to reconstruct 



the parton kinematics in heavy ion collisions as well and, t 



partonic energy loss [k]]. Recent phenomenological analysis 



rerefore, to establish the pattern of 
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13| indicate that studying distri- 



butions in the fragmentation fraction variable z is particularly sensitive to the medium modified 



3 



fragmentation functions at RHIC. 

Furthermore, a similar observable corresponding to the production of a jet and a photon can 
be used, as proposed in [I5I ]. to constrain the photon fragmentation functions at RHIC in a region 
barely accessible by LEP kinematics. 

In order to make reliable quantitative predictions for a high-energy process, it is crucial to 
determine the NLO QCD corrections to the Born approximation. In hadronic collisions, cross 
sections computed at the lowest order in perturbation theory are only of qualitative value, since 
higher order corrections are usually sizable, in some cases as large as or even larger than the Born 
contribution. Besides that, the LO calculation is severally affected by the dependence on the 
'"unphysical"' factorization and renormalization scales, dependence that can be partially cured 
only by including the higher order corrections. Also, for less inclusive processes, the inclusion of 
extra partons in the NLO perturbative calculation allows one to improve the matching between 
the theoretical result and the realistic experimental conditions. Furthermore, the usual '"naive"' 
relations between the momentum of the outgoing particles and the kinematics of the participating 
partons in the hard process are usually exact only at leading order accuracy, being mandatory 
to establish their range of validity when higher order corrections are taken into account. NLO 
corrections for 7 + h were computed for unpolarized collisions in [jjj], where the corrections are 
found to be large, on the order of 70%, even for LHC energies. 

The purpose of this paper is, therefore, to extend the available calculations by computing the 
NLO corrections to hadron plus prompt-photon production in polarized hadronic collisions and 
to analyze its phenomeno logical impact. Nowadays there are several versions of the subtraction 



method that allow the computation of any infrared-safe quantity to NLO accuracy. We re 



the formalism of Ref. 



161 ]. which has been used previously to compute observables in jet 
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yon 



prompt-photon [19], single-hadron [20J and hadron+jet production 2l|], for both polarized and 
unpolarized cases. 

Most of the advantages for detecting a photon discussed above rely on the fact that the photon 
couples directly (pointlike) to the partons in the hard process and propagates as a colorless state. 
Unfortunately, this is true only for direct photons, and the cleanliness of the process is spoiled 
by resolved photons that can be produced from a parton in the QCD reaction that fragments in 
a photon plus a number of hadrons. Several isolation algorithms have been proposed to reduce 
this contribution. The main idea in all of them is that nondirect photons are surrounded by 
hadronic debris, so measuring the amount of hadronic energy around a photon allows us to define 
a selection criteria. In this paper we apply the isolation prescription introduced by Frixione in 
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Ref. 22], which has the theoretical advantage that, still being infrared-safe, completely eliminates 
the unwanted resolved component to the cross section. 

Relying on this isolation algorithm, we construct a Monte Carlo code that calculates any 
infrared-safe observable to NLO precision in hadron plus direct photon production. The code 
provides results for both unpolarized and polarized collisions and allows us to impose experimental 
cuts and to perform a detailed phenomenological analysis. Specifically, we use the Monte Carlo 
code to explore the sensitivity of some observables on the spin-dependent gluon distributions in 
polarized proton-proton collisions. Also, we study correlations between certain sets of external 
experimentally-accessible variables and real partonic momentum fractions in the unpolarized case. 

The paper is organized as follows: in Sec. II we describe the theoretical tools used to perform 
the NLO calculation and study the perturbative stability of the LO and NLO results. In Sec. Ill 
we analyze some general aspects of unpolarized and polarized collisions. In Sec. IV, we focus on 
the polarized case, showing predictions for asymmetries at different center-of-mass energies and 
both RHIC experiments. Finally, we present the conclusions in Sec. V. 



II. NLO CORRECTIONS AND PERTURBATIVE STABILITY 

The factorization theorem allows one to separate the nonperturbative effects associated 
with the hadronic structure from those calculable within perturbation theory and, therefore, to 
write the hadronic cross section as 

da^ H '^\K 1 ,K 2 ,K 3 ,K A )= J dx 1 dx 2 dzfl*\x l , f i J )fg>\x 2 ,» I )D$Xz, t i F ) 

da aia2a3J (xiKx, x 2 K 2 , K 3 /z, K 4 ; fi R , /x 7 , /i F ) 

+ J dx 1 dx 2 dz 1 dz 2 f^\x 1 , t i I )f^\x 2 ,fi I )D^(z 1 ,fi F )D^(z 2 , t i F ) 

0,10,20,30,4 

d^a 1 « 3 «3o 4 (iCi^i, x 2 K 2 ,K 3 /z 1 , K 4 /z 2 ; fx R ,fXj,n F ) , (1) 

where H\ and H 2 are the colliding hadrons, K\ and K 2 their four-momentum, and the sum is 
performed over the parton flavors which can contribute to the process. Here, K3 is the final state 
hadron four-momentum, while K4 denotes the photon momentum. The first term represents the 
direct contribution whereas the second one accounts for the resolved component Q. da aia2a31 and 
do~ aia2a3a4 are the corresponding subtracted partonic cross sections, in which collinear singularities 
1 Notice that in general the decomposition between direct and resolved contributions from NLO accuracy is not 
physical. Nevertheless, with the isolation prescription used along this work only the first term survives, which 
implies that the direct contribution to the total cross section can be unambiguously defined. 
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have been removed through the introduction of suitable counterterms. fL is the partonic distri- 
bution function (PDF) for a parton cij inside an hadron Hi, while D^* denotes the fragmentation 
function associated with a parton a$ which hadronises into a hadron h. In the resolved term also 
the necessary Dap photon fragmentation function appears. The expression in Eq. (pQ) implicitly 
assumes that the photon and the final state hadron are well separated in angle, to avoid extra 
collinear configurations that would require the introduction of double fragmentation functions 
24}]. Here, /// and //p correspond to the initial and final state factorization scales, while (Ir is 
the renormalization scale. The full set of NLO QCD corrections to both direct and resolved con- 
tributions were first computed for unpolarized collisions in Ref. [jjj]. The analogous of §B) for 
polarized cross section is obtained by replacing the parton distributions and the partonic cross 
section by its polarized expressions, A/^. and dAa, respectively. As usual, the longitudinally 
polarized asymmetry is defined by the ratio between the polarized and unpolarized cross sections 

, h dAa . . 

^ = -s- P) 

As stated in the introduction, resolved contributions take into account those events in which 
the photon is generated from parton fragmentation. In order to eliminate its contribution, we rely 
on the photon isolation prescription introduced in Ref. (22]. This algorithm can be summarized 
in four steps: 

• Identify the photon in the final state and define a cone surrounding it with a radius r in 
the rapidity-azimuthal plane; 

• If all QCD partons lie outside the cone, then the photon is isolated; 

• If there are QCD partons inside the cone, calculate their distance to the photon 



which is the total transverse hadronic energy accumulated in a cone of radius r around the 
photon; 

• Introduce an arbitrary function £(r) satisfying lim r ^o£( r ) — 0. Then the photon will be 
isolated if Ex{r) < £(?"), for every r < tq. 




(3) 



and define 




(4) 
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Notice that this prescription does not forbid the emission of soft gluons in any region of the phase 
space, but it does eliminate the divergences arising from the kinematical configuration where 
a quark becomes collinear to the photon. For this reason, the prescription does not spoil the 
cancellation of infrared singularities and can be applied to any infrared-safe observable. 

After the implementation of the isolation criteria only the direct contribution to the cross 
section survives as 

da^ H ^ h ^ ISO \K 1 ,K 2 ,K 3 ,K 4 )= J dx 1 dx 2 dzfW(x 1 ,» I )fg>\x 2 ,» I )DW(z,ii F ) 

a l a 2 a 3 

d^^ixiKx, x 2 K 2 , K 3 / z, K 4 ; fi Rl fij,fi F ) , (5) 

where the (ISO) superscript indicates that the isolation algorithm has been included in the cor- 
responding measurement function. Following the procedure detailed in Ref. jit] ], the calculation 
is implemented in a Monte Carlo-like code that allows us to compute any infrared-safe observable 
for the process HiH 2 — > hj. It is worth noticing that the same code computes both unpolarized 
and polarized cross sections, since the formal structure of the corrections is exactly the same. 

In this work, we mostly concentrate on the phenomenology of charged pion production ac- 
companied by a back-to-back photon. Our default results will correspond to the kinematics of 
the PHENIX experiment at RHIC with a center-of-mass energy of i/i = 200 GeV, but we will 
also present asymmetries at \fs = 500 GeV and for the kinematics corresponding to the STAR 
experiment as well. Unless otherwise stated, we require the pion transverse momenta to be larger 
than 2 GeV and the one for the photon to be restricted to the 5 GeV < ]P T < 15 GeV range. These 
cuts are imposed to take into account experimental limitations to identify direct photons [?], sj]. 

The rapidities of the pion and the photon are limited to the range \r]\ < 0.35, according to 
the constrains imposed by the PHENIX detector q The particles have to be separated in the 
azimuthal angle by A0 = I0 71 " — <p\ > 2 to ensure the pion is produced from a parton in the 
'opposite side' to the photon. For the isolation algorithm we chose the radius as r = 0.4, the 
functional form 



1 — cos(r ) 

and set the parameters in Eq. by n = 4 and e 7 = 1. 

Since there are two different hard scales in the process, given by the transverse momenta of 
the pion and the photon, the default initial (fij) and final (fip ) state factorization scales and the 



2 For the sake of simplicity we assume full azimuthal coverage for our phenomenological results. The rapidity 
limits are set to \r/\ < 1 for STAR. 
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renormalization scale (//r) are chosen as their average. Therefore, unless otherwise stated the 
scales are defined as 

_ Pt + Pt />t\ 

Hi = Hf = Hr = H = . (7) 

In the unpolarized case, we rely on the MRST2002 parton distribution functions (PDFs) 251 ]. 
since it is the reference set for the DSSV polarized PDFs that we will use as the default density 
in the computation of the corresponding asymmetry. For the fragmentation functions we rely on 



the DSS 



26] set that provides full flavor and charge separation at NLO. We will later discuss the 



consequences of using different sets of parton distributions and fragmentation functions. 




FIG. 1. Unpolarized (solid line) and polarized (dashed line) NLO K-factors for the pion transverse 
momentum distribution (left) and the variable Ze (right). The choice of the factorization and renormal- 
ization scales corresponds to the default value [ij = fjLp = fJ-R = fi= (p T +Pt) The P°l ar "i ze d K-factor 
is not presented in the right-side plot since the polarized cross section has nodes in that range of Ze- 

It is important to remark that we will always use NLO PDFs and fragmentation functions 
for both LO and NLO calculations. This choice is particularly important in the polarized case 
where the uncertainty in the distributions can render artificially large corrections because of the 
differences between the LO and NLO polarized densities. In other words, any difference we find 
between LO and NLO cross sections will arise only from the corresponding partonic cross section. 
Also, we use in both cases the two- loop expression for a s , and the one- loop correction to a. 
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A useful way to quantify the size of the QCD corrections is by means of the .fT-factor, defined 
as the ratio between the NLO and LO results. In Fig. [T]we show the i^-factor as a function oi p T 
and of the scaling variable Ze [] defined as 

Z E = - cos (8) 

Pt 

where (ff and (jf are the azimuthal angle associated to the pion and the photon, respectively. We 
will discuss in the next section the usefulness of this variable, which at lowest order represents 
the argument (z) of the fragmentation function in Eq. (jSJ). As it happens in the case of single- 



0.01 



1E-3 



1E-4 



CD 

-Q 1E-5 



1E-6 



1E-7 



1E-8 



NLO 



LO (x0.1) 



4=2 

•5=1 
4=0.5 



NLO 



: LO(x0.1) 



- E 



CM 200 GeV, h|<0.35 



FF DSS n 

PDF MRST2002 NLO 
5GeV<p! r <15GeV 



p" T (GeV) 



0.01 



- 1E-3 



1E-4 



1E-5 



1E-6 



'I I I I I I I I I I I I I I ' 

34567893456789 



-. 1E-7 



1E-8 



FIG. 2. Scale dependence of the unpolarized cross section in terms of the transverse momentum of the 
pion. The left side corresponds to the simultaneous variation of factorization and renormalization scales 
as \ii = fj,p = [xr = £fj, = £ (pj, + Pt) /% while the right side shows the result when the renormalization 
scale is kept fixed to the default value. 



hadron production [20, |27j, the NLO corrections turn out to be very large, with a i^-factor that 
lies in the range 2.2 < K < 3.2, and with more moderate (but still sizable) corrections for the 
polarized cross section at lower transverse momentum. Notice that for these observables the K- 
factor is a nontrivial function of p T and Ze, respectively. Furthermore, the pattern for the QCD 



This variable is denoted as Xb in Refs 
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corrections in the polarized case considerably differs from the unpolarized one, which will result 
in non-negligible corrections for the asymmetries. 

Besides the if-factor study, we need to estimate the corresponding theoretical uncertainty for 
this calculation. Because we are not able to evaluate the full-perturbative expansion, we can 
indirectly test the quality of the NLO expansion through a scale-dependence analysis. As it 
is well known, factorization and renormalization scales are not physical parameters of the real 
problem, so the exact result must not depend on them. But theoretical predictions do have such 
a dependence, arising from the truncation of the perturbative expansion at a fixed order in the 
coupling constant a s . A large dependence on the scales, therefore, implies a large theoretical 
uncertainty. A full NLO analysis of scale dependence for unpolarized collisions was presented in 
Ref. 141 ] . Nevertheless, and for the sake of completeness, before showing the new results for the 
polarized cross section we recall the main features in the unpolarized case at RHIC kinematics. 
As a first approach to analyze the perturbative stability of the observable through the scale 
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FIG. 3. Same as Fig. [H but for the polarized cross section. The plot starts in this case at pj, = 6 GeV 
since at lower values the polarized cross section changes sign. 

dependence, we set again all factorization and renormalization scales to be equal and vary them 
by a factor of 2 up and down with respect to the default choice, i.e., //j = \ip = Hr = £ {j>t + Pr) 
with £ = 1/2, 1, 2. The left side of Fig. [2] plots the corresponding scale dependence of the LO and 
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NLO unpolarized cross sections. When the three scales are varied together, both LO and NLO 
results exhibit a very similar 0(±5O%) variation, which is at first sight disappointing since one 
expects an improvement after including the NLO corrections. However, this is not an indication 
of a failure of the perturbative expansion as it might seem. On one hand, a large scale de pen dence 
is not truly unexpected since similar results are found for single-hadron production 20j,|27j|. The 



fact that no visible improvement is found at NLO for this choice of scales has also been observed 
in the case of inclusive prompt-photon production 19|, which relies on the same matrix elements 
used for this calculation. This trend can be attributed to the particular choice of scales, produced 
because of an incidental cancellation between the effects of the factorization and renormalization 
scale variations at LO when they are varied simultaneously. 

Therefore, a better way to analyze the issue of perturbative stability in this case is to perform 
an independent variation of the scales. As an example, we keep the renormalization scale fixed to 
the default value (i.e., li r = fi) and vary only the initial and final state factorization scales in the 
way described above (/!/ = llf = £a*)- The result of this study is shown in the right side of Fig. [2} 
where we do observe that the LO cross section presents an average fluctuation of ± 30% whereas 
the corresponding NLO result varies only about ± 15%. Thus, in this scenario we can assert that 
NLO is truly improving the accuracy of the calculation. The dependence on the renormalization 
scale, when the factorization scales are fixed to the default value, remains almost the same at 
LO and NLO. The lack of reduction on the lir dependence at NLO is due to the isolation cuts 
which, while being an infrared-safe prescription, still perturb the cancellation of the soft-gluon 



effects and therefore have an impact on the renormalization scale dependence [28|. Most of the 
scale dependence, particularly the one arising from the renormalization scale, will cancel in the 
cross section ratios built in order to analyze the nuclear effects, and in the asymmetries defined 
in polarized collisions as well. Nevertheless, despite the improvements at NLO that confirm the 
perturbative stability, it is clear that the theoretical uncertainty remains rather sizable for the 
cross sections themselves Q. All the results in Fig. [2] correspond to the pion transverse momentum 
distributions. We observe similar features when we study the scale dependence of the unpolarized 
cross section as a function of the scaling variable Ze- The trend of the scale dependence turns 
out to be the same when the variations are performed around a different default choice, like the 
invariant mass of the photon-pion pair or different combinations of their transverse momentum 
like pj-, pf. and max(pj,pj). In particular, we do not find any reasonable choice that presents a 
4 Variations of the order of 40% are quoted in Ref. [14[ for the ir-f invariant mass distributions at the LHC after 
implementing a different isolation prescription. 
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visible decreases in the scale dependence. 



1E-3 



1E-4 



03 



1E-5 



1E-6 



Unpolarized 

PDF MRST2002 NLO 



E CM =200 GeV 

FF DSS n 

p;>2 GeV, h|<0.35 

5GeV<p!<15GeV 



1 ■ i 1 ■ 

0.03 0.06 



0.09 



0.12 



0.15 



T 
0.18 0.03 



X' 



Polarized 

PDF DSSV NLO 



■4.0x10 



-2.0x10 



; -0.0 



0.06 



0.09 



0.12 



0.15 



0.1E 



> 
q 



cr 



-2.0x10 3 



-4.0x10 



-6.0x1 0" 7 



FIG. 4. NLO cross section as a function of X^ ca \ (solid lines), Xr cc (dashed lines) and X^ cc (dotted 
lines), for the polarized (right) and unpolarized (left) cases. Factorization and renormalization scales are 
set to the default value. 



For polarized collisions, the reduction in the scale dependence at NLO is much more evident, 
as shown in Fig. |3] This feature has been observed for many observables in hadronic collisions 
and can be mainly attributed to the less singular behavior of polarized parton distributions and 
the corresponding evolution kernels that drive the initial state factorization scale-dependence. 



III. PHENOMENOLOGY 

As discussed in the introduction, counting with the photon and pion kinematics allows us to 
attempt for a reconstruction, to some degree of accuracy, of the partonic momentum fractions x\$ 
and z. It is worth noticing that at the Born level, the momentum fractions X\i can be written as 

f T exp {rf) + pi exp (r/ 7 ) 

-<M,Rec — "F l y J 

V s 

_ p£exp(-7f) +p£exp(-rf) 

^-2,Rec — 7= , l-LUJ 
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or equivalently at this level, as defined by the experimental collaborations, 

p], exp (jf) — cos (4> n — 7 )pj exp {rp) 



x' 



Rec v/i 

]P T exp (— r/ 71 ") — cos (0 71 " — (p)p^p exp {—f] 1 ) 



2,Rcc 



(11) 
(12) 



At the Born level, these variables are equal to the initial partons momentum fractions, x% and x^ 
respectively, because only the 2—^2 partonic process enter in the calculation and cos (0 71 " — 7 ) = 
— 1 for the trivial back-to-back configuration. However, NLO computations require the inclusion 
of virtual and real corrections, the last being associated with 2—^3 partonic processes. Therefore 
Xj R ec and X- Rec may differ from the real momentum fraction Xi and between themselves. Counting 
with the Monte Carlo code presented in the previous sections, we can analyze to what extent these 
naive relations are modified by the NLO corrections. Recalling the invariance of the physical 
problem under the exchange of partons 1 and 2 at RHIC, we study the cross section distribution 
in terms of the reconstructed variables Xr gc = X^Rec, X' Kec = X[ Rec and also the true momentum 
fraction X Real = x x . 
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FIG. 5. Correlation for Xnec (left) and X^ cc (right) versus X^^. We use a gray scale picture to indicate 
the probability of finding an event with a certain value of X Roa j and X^ ec (or X^ cc ): the darker the bin, 
the larger the correlation between the variables. For the sake of presentation we plot values of X multiplied 
by 1000. 



As we can observe in Fig. HJ the shape of the NLO cross section as a function of these 
three variables is rather similar for both the unpolarized and polarized cases in the kinematical 
region, where the corresponding cross section is sizable. For the unpolarized result, da/dX^ ec 
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seems to be closer to da/dX^ eSb \ than da/dX-& ec ] in the polarized case, using the DSSV NLO 
parton distributions, the trend is the opposite. We also observe from Fig. @] that, with the 
kinematical cuts introduced above, the dominant range in x for this observable is determined by 
0.05 < x < 0.15. 
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FIG. 6. "'Gluonic ratio'" R g as a function o/ Xr cc (left), p^> (center) and Ze (right), for unpolarized 
proton-proton collisions at RHIC. We analyze separately the cross sections for positive (solid lines) and 
negative (dashed lines) pions in the final state. 



In any case, comparing the shape of the cross section as a function of different variables is 
not enough to ensure that these variables have a strong correlation. Since the Monte Carlo 
implementation of the NLO calculation provides information about the kinematics of all the 
particles involved in the process, we can perform an event-by-event comparison of the values of 
X Rec , X Rec and X Real , weighting each possible configuration with its associated cross section. The 
result is presented in Fig. El where a gray-colored scale is used to indicate the probability of 
having an event with such values of the variables. A perfect correlation would be described by a 
black diagonal, as it occurs to LO accuracy: as we can observe X Rec and X Rec are both strongly 
correlated with X Rea i. Using the definition of correlation for a nonuniform probability space, we 
find that the correlation coefficient is larger than 0.76 for both pairs of variables (X Rea i, Xr 6C ) and 
(XR ea i, X Rec ), when the bin size is set to 0.016. 
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FIG. 7. Contributions from the different initial state partonic channels to the polarized cross section 
as a function of pj, (left) and X^ec (right), for polarized proton-proton collisions at RHIC. Here qQ 
corresponds to the sum of all channels initiated only by quarks or antiquarks. 

The reconstruction of the initial state momentum fractions might become very useful in order 
to extract information about parton distributions, and specially the polarized gluon density in 
polarized collisions. Regarding that, it is interesting to study the dependence of the cross section 
on the different partonic channels, to quantify the sensitivity on the gluon distribution. For that 
purpose we define the '"gluonic ratio"' as 

dagg + da qg 
tig = — -r , 

"OTotal 

which measures the relative importance of partonic processes initiated by at least one gluon. 

In Fig. [6] we plot R g as a function of pj, Ze and Xr 6C , where we distinguish between positive 
(solid lines) and negative (dashed lines) pions. Positively charged u quarks and gluons dominate 
the inner structure of the proton, so it is more probable to have gluons or u quarks in the final 
partonic state, which involves a higher probability of a fragmentation into positive pions and 
increases the sensitivity on the gluon distribution for tc + + 7 compared to tt~ + 7. As can be 
observed in Fig. [6l R g (iT + ) exceeds 0.80 for a wide kinematical range. Particularly, in the range 
of 3 GeV < Pt < 10 GeV, more than 88 % of the events contain at least one gluon in the 
initial state, confirming that hadron + photon production in hadronic collisions might provide a 
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clear source of information on the gluon distribution. Of course, these ratios will considerably 
change in the polarized case due to the modification in the partonic cross sections as well as the 
different behavior and magnitude of the polarized distributions. As an example, we show in Fig. 
[7] the contribution from each initial state partonic channel to the polarized cross section, for the 
(positively charged) pion transverse momentum and X^ ec distributions. It can be observed that 
the gluon contribution, mostly due to the qg channel, is rather sizable in the polarized case, even 
when a set with small gluon polarization like the DSSV set is used. In the next section we will 
analyze the sensitivity on Ag in more detail. 
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FIG. 8. NLO cross section as a function of Z^ cai \ (solid), Zr oc (dash) and Ze (dot), for the unpolarized 
(left) and polarized (right) cases. Factorization and renormalization scales are set to the default value. 

On the other hand, in the case of heavy ion collisions it is of special interest to get a grip on the 
fragmentation variable, which would allow one to estimate the initial energy of the parton that 
undergoes hadronization and, eventually, loses energy when traveling through the high density 
matter formed in the collision. 

In theprevious section we already introduced the variable Ze (see Eq. (jHJ)), which was employed 
in Ref. [7| to study partonic energy loss in Au-Au collisions at RHIC. At the Born level, Ze is 
equal to the true momentum fraction of the pion originated from the hadronization of parton 3, 
which we denote as -ZR ea i = z. Alternatively, one might define 

ZRcc = ^ > ( 14 ) 

Pt 
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which is also equal to Ze at LO. As happened with X^ ec and A Rcc , Z^ ec and Ze can be used 
to approximate the real momentum fraction Z^ ea \ m a NLO computation. To test the validity 
of this estimation, we plot in Fig. Elthe cross section as a function of ZR ea i (solid), Zr cc (dash) 
and Z E (dot), for the unpolarized (left) and polarized (right) cases. The three functions have a 
similar shape, although the discrepancies become larger for the polarized case in the region 0.15 
< z < 0.45, where the polarized cross section changes sign when the DSSV set of polarized PDFs 
is used. It is worth noticing that the cross section is dominated by the fragmentation of a parton 
into a pion with a momentum fraction in the range 0.1 < z < 0.6. 
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FIG. 9. Correlation for Z^ ec (left) and Ze (right) versus ^R ea i- As in Fig. [5|, we use a gray scale picture 
to indicate the probability of finding an event with a certain value of ivReal and Zr cc (or Ze): the darker 
the bin, the larger the correlation between the variables. 

As we did for the initial state momentum fractions, we also performed an event-by-event study 
to quantify the correlation of the pairs of variables (Z Real , Z Rec ) and (Z Re3j i,Z E )- Applying the 
same algorithm used for initial state momentum fraction, the correlation coefficient is 0.77 for 
(ZReai, Zrcc) and 0.71 for (ZR ea i, Ze), when the bin size is set to 0.08. Therefore, we conclude 
that Z^ cc allows for a more accurate reconstruction of ZR ea i at NLO accuracy. In spite of having 
a high correlation coefficient, the reconstruction is not exact and both approximations show a 
bifurcation in their trends as z increases, due to the fall of the cross section. This situation can 
be observed in Fig. [91 where we can appreciate that the deviations from the exact values are 
anyway typically smaller than 15%. These correlations can be improved by setting stronger cuts 
on the transverse momentum of the photon, like > 10 GeV, but with the corresponding loss 
in cross section. This level of reconstruction might help to obtain a better understanding of the 
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energy loss of partons when traveling through the dense medium created in heavy ion collisions 
and, particularly, to test whether the parton energy loss can be accounted for via an effective 
modification of the fragmentation functions. We have analyzed other possible definitions for the 
reconstructed version of the variables x and z, all of them equivalent to LO accuracy. We do 
not find sensible differences in the corresponding correlations with the real variables compared to 
those already presented here. 
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side plot) charged pion computed using different set of parton distributions: MRST2002 (solid lines), 
CTEQ6M (dotted lines), CT10 (dashed lines) and MSTW2008 (dotted- dashed lines). 



Before moving to the phenomenological analysis of the asymmetries, we would like to quantify 
the possible uncertainties on the analyzed observables due to the unpolarized parton distributions 
and fragmentation functions. As stated before, we rely on MRST2002 as the default set because 
the DSSV analysis is based on those distributions. In Fig. [TO] we show the X Rec distribution for a 
photon accompanied by a positively (left-side plot) and negatively (right-side plot) char ged pion 
computed using different sets of parton distributions: MRST2002, CTEQ6M CT10 |3o| and 



MSTW2008 



31| . While the differences become more visible in the kinematical range where the 



cross section becomes rather small, in general we do not observe variations large enough to spoil 
the interpretation of the measured asymmetries. This is not unexpected, since modern sets of 
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FIG. 11. Ze distribution for a photon accompanied by a positively (left-side plot) and negatively (right- 
side plot) charged pion computed using different set of fragmentation functions: DSS (solid lines), AKK 
(dashed lines) and Kretzer (dotted lines). 



PDFs provide distributions that agree within a few percent between them. The situation is more 
complicated when one turns to the fragmentation functions, since those distributions still suffer 
from much larger uncertainties. In order to study the dependence on the set of fragmentation 
functions used, we plot in Fig. [TT] the Ze distribution for both positively and negatively charged 
pions obtained using the DSS 26], AKK 32| and Kretzer {33]] sets of fragmentation functions. 
As anticipated, the variations observed here are more sizable than the ones found for the parton 
distributions. Nevertheless, we will demonstrate in the next section that the dependence on the 
fragmentation functions, which contribute in a similar way to both polarized and unpolarized cross 
sections, partially cancels when computing the corresponding asymmetries. 



IV. SPIN ASYMMETRIES 

As a way to directly analyze the sensitivity of the process on the polarized gluon distribution, 
we will compute the NLO asymmetries with 3 different sets of spin-dependent densities: DSSV 
GRSV (standard) (3^ and GS-C [33|. The corresponding NLO distributions at Q 2 = 50 GeV 2 , 
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FIG. 12. Polarized gluon density at Q 2 = 50 GeV 2 from different sets of polarized PDFs (left) and their 
ratios to the unpolarized distribution (right). 



a typical scale for this process, are shown in the left-side of Fig. [121 As can be observed, the 
expectations from the three sets are quite different. Here, the DSSV distribution corresponds to 
the best fit from the latest global analysis of all polarized data [1], while the GRSV set can be 
considered as an (overestimate of the) '"upper bound"' for the allowed range of gluon densities, 
even though this set has already been ruled out by the available data. The GS-C set provides a 
distribution compatible with the requirement of a small gluon polarization in the range 0.05 < 
x < 0.3 but with a node in that region and a very different behavior at smaller x compared to 
the DSSV set. The right side in Fig. [12] shows the ratio between the corresponding polarized 
distribution and the unpolarized MRST2002 set. In the quark sector, the dominant distributions 
are very similar among the different sets. Therefore, we can expect that any differences between 
predictions for the asymmetries that are found when using different polarized parton density sets 
are to be attributed to the sensitivity of the observable to Ag. 

For the presentation of our phenomenological results we will concentrate on the production of 
positively charged pions, since that observable provides, in the polarized case, a stronger sensitivity 
on the polarized gluon distribution. We start with the expected asymmetries by looking first at 
7r± + 7 production in terms of the transverse momentum of the pion as shown in Fig. [13j As 
expected, the asymmetries for the DSSV and GS-C distributions turn out to be rather small. 
When a set with a larger gluon distribution is considered, like GRSV, the asymmetries increase 
to the ~ 2% level for 7r + production. A rough estimate of the minimum value of the asymmetry 
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FIG. 13. Expected NLO and LO asymmetries for n + (left) and ir~ (right) production at RHIC in terms 
of the transverse momentum of the pion for different sets of polarized PDFs: DSSV, GS-C and GRSV. 
The histogram in the left- side plot represents the value A m i n as defined in Eq. U5\) . 

observable at RHIC due to statistics can be obtained by the well-known formula 

Amm = ^^hre' (15) 

where £ is the integrated luminosity, P is the polarization of the beam, and the factor e < 1 
accounts for experimental efficiencies; a is the unpolarized cross section integrated over a small 
range in transverse momentum of the pion. The quantity defined in Eq. (|15p is plotted (histogram) 
in the left side of Fig. [131 for e = 1, P = 0.7, C = 100 pb -1 , and a pr-bin size of 2 GeV. From the 
value obtained with these particular set of parameters one can notice that the minimally observable 
asymmetry at RHIC is of the order of magnitude of the expected value, and typically larger than 
the one obtained using the DSSV set. A similar situation has been observed for the asymmetries 
already measured at RHIC in the case of jet or single-inclusive hadron production, which turned 
out to be the key measurements to determine the polarized gluon distribution so far [l|. Therefore, 
one might expect that the measurement of ir + + 7 at RHIC will also play an important role in 
future global analysis. As anticipated in the unpolarized analysis (see Fig. [7]), the sensitivity on 
the gluon distribution drops considerably for tt~ , an effect which is reflected in the decrease of the 
corresponding asymmetry. Nevertheless, the separate measurement of asymmetries for positively 
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and negatively charged pions can provide extra information on the polarized quark flavor structure 
once incorporated in a global analysis of all observables. It is also worth noticing from Fig. [TS] 
that the NLO corrections in general tend to reduce the asymmetries, but the effect is nontrivial 
and strongly depends on the transverse momentum of the pion and the polarized set of parton 
distributions used. 
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FIG. 14. Same as Fig. [73] but in terms of the scaling variable X-R ec (left) and Ze (right). 



A similar study is presented in Fig. HU where we analyze the corresponding asymmetries for 
7T + + 7 production in terms of the scaling variables Xr S c (left) and Ze (right). In particular, the 
left-side plot in Fig. [TH reflects the x-shape and order of the curves for the AG/G ratio plotted 
in the right side of Fig. [121 confirming the strong sensitivity of this observable on the gluon 
polarization. 

As discussed in the previous section, one issue that can veil the possible extraction of infor- 
mation about the polarized parton distributions arises from the uncertainty in the fragmentation 
functions. A comparison at the level of the asymmetries in terms of both scaling variables X^ ec 
and Ze is presented in Fig. [151 where we confront the results obtained using the DSS (left-side 
plots) and AKK2008 (right-side plots) of fragmentation functions to NLO accuracy. For most 
of the kinematical range that can be covered by the RHIC experiment, the dependence in the 
fragmentation functions, quite visible at the level of the cross section in Fig. [HI almost cancels in 
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FIG. 15. NLO asymmetries in terms of the scaling variable Ze (up) and X-^ ec (down). The left-side 
plots correspond to results obtained using the default DSS set of fragmentation functions while those in 
the right were obtained using the AKK2008 set. 

the asymmetry ratios, with the exception of the end points in Ze, where fragmentation functions 
are hardly constrained by the available data. 

The results presented so far correspond to a center-of-mass energy y/s = 200 GeV and were 
restricted to the kinematics accesible to the PHENIX experiment. In the following we discuss 
the extension of the phenomenological analysis of the asymmetries to a center-of-mass energy 
y/s = 500 GeV and also account for possible measurements at STAR, which extends the rapidity 
coverage to \t)\ < 1. 

In Fig. HE] we present the asymmetry for 7r + production at y/s = 500 GeV in terms of the 
transverse momentum of the pion for PHENIX (leftside plot) and STAR (rightside plot). As 
can be observed, the asymmetries show the same pattern found for y/s = 200 GeV but almost a 
factor of 2 smaller in magnitude. For the same value of transverse momentum, the cross section 
becomes sensitive to smaller values of x when the center-of-mass energy increases. Since the 
polarized parton distributions fall faster than the unpolarized ones at small x, that results in the 
reduction observed in the asymmetries. Also, at smaller x the quark-initiated channels play a 
more important role in the polarized cross section and tend to reduce the gluon sensitivity of the 
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FIG. 16. Expected asymmetries for tt + production at RHIC in terms of the transverse momentum of the 
pion for center- of-mass energy \fs = 500 GeV and PHENIX (left side) and STAR (right side) kinematics. 
The histogram represents the value A m i n as defined in Eq. $15\). 



asymmetry. 

For this particular observable one can not appreciate large deviation between the predictions 
for PHENIX and STAR. The integration over a larger range in rapidity at STAR provides a 
greater cross section, with the corresponding increase in statistics, but that normalization effect 
almost cancels in the asymmetry. We also include in Fig. [16] the estimate of the minimum value 
of the asymmetry observable at RHIC for both experiments obtained using Eq. (fT5l) . assum- 
ing an integrated luminosity of C = 300 pb -1 . We observe that, despite the decrease in the 
asymmetries because of the larger center-of-mass energy, the analyzing power of the observable 
actually improves due to the much higher statistics that can be collected at this energy for both 
experiments. 

Finally, we present in Figs. [T7] and [18] the asymmetries for tt + production in terms of the 
scaling variables Xr 6C (left-side plots) and Zg (right-side plots) for PHENIX and STAR kinematics, 
respectively. In the particular case of the X^ ec distribution, we observe that one can span values 
of x about a factor of 2 smaller than at \fs = 200 GeV. Furthermore, one can appreciate subtle 
differences between the results for PHENIX and STAR, mostly due to the fact that the larger 
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FIG. 17. Same as Fig. [77] for center- of-mass energy y/s = 500 GeV and PHENIX kinematics. 

rapidity coverage of STAR slightly broadens the range in the momentum fraction of the initial 
state partons that contribute to the observable. 



V. CONCLUSIONS 



In this paper we presented the NLO corrections for the unpolarized and polarized cross sections 
for the production of a hadron and a back-to-back isolated photon. 

The corrections are found to be nontrivial: '/\~-factors' are larger for the unpolarized cross 
section than for the polarized one, typically resulting in a reduction of the asymmetry at NLO. 
It is shown that the perturbative stability of the hadron plus photon cross section improves after 
including the NLO contributions, particularly when analyzing the dependence on the factorization 
scales, but the total scale dependence remains rather sizable due to the implementation of isolation 
cuts for the photon. 

We confirmed that '"naive"' LO relations between the photon and pion momentum retain a 
strong correlation to the true partonic fractions which are the arguments of the parton distributions 
(xi and X2) and fragmentation functions (z) when NLO effects are taken into account. Counting 
with more precise distributions on the fragmentation fraction z might permit one to obtain a better 
understanding of the energy loss of partons when traveling through the dense medium created in 
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FIG. 18. Same as Fig. ^T^for STAR kinematics. 



heavy ion collisions. 

Finally, the possibility of looking at a charged pion accompanied by a back-to-back isolated 
photon is studied phenomenologically in detail for polarized collisions, as feasible at RHIC. We 
find that the asymmetries for n + production, in terms of both dimensionless variables X^ ec and 
Ze-i are sensitive to the polarized gluon density in the range 0.05 < x < 0.15 and might contribute 
to a better understanding of the spin content of the proton. 
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